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Abstract 

The fluorescence yield for dry air and pure nitrogen excited by electrons is calcu- 
lated using a combination of well-established molecular properties and experimental 
data of the involved cross sections. Particular attention has been paid to the role 
OO I of secondary electrons from ionization processes. At high pressure and high energy, 

Oj ' observed fluorescence turns out to be proportional to the ionization cross section 

.^ . which follows the Born-Bethe law. Predictions on fluorescence yields in a very wide 

^^ I interval of electron energies (eV - GeV) and pressures (1 and 1013 hPa) as expected 



from laboratory measurements are presented. Experimental results at energies over 
1 MeV are in very good agreement with our calculations for pure nitrogen while 
discrepancies of about 20% are found for dry air, very likely associated to uncer- 
tainties in the available data on quenching cross sections. The relationship between 



^ ■ fluorescence emission, stopping power and deposited energy is discussed. 
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1 Introduction 



The fluorescence radiation from air showers generated by ultra-high energy 
cosmic rays UHECRs provides a precise determination of the longitudinal 
profile allowing the reconstruction of the primary properties. This technique [1] 
was first successfully used by the Fly's Eye experiment [2] and later by HiRes 
[3]. Fluorescence telescopes are being used by the Pierre Auger Observatory [4] 
to operate simultaneously with a giant air shower array. The planned Telescope 
Array project [5] also relies in this technique. On the other hand, the satellite- 
based experiments EUSO [6] and OWL [7] are being designed for the detection 
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of the fluorescence traces of air showers viewing downward from the top of the 
atmosphere. 



The HiRes collaboration has reported measurements of the energy spectrum of 
UHECRs [8] which are in disagreement with those of the AGASA air shower 
array [9]. In order to achieve a reliable energy calibration of a fluorescence 
telescope, accurate values of the air fluorescence yield are required. This need 
has promoted a number of experiments for the measurement of this parameter 
in the wavelength interval of interest in this technique, i.e. 300 - 400 nm. As 
is well known, in this spectral interval, the fluorescence light comes from the 
First Negative System of N^ (B^S+ -^ X'^T,^) and the Second Positive System 
of N2 {C^Hu -^ B^IIg) which in the following will be called IN and 2P systems 
respectively. 



Kakimoto et al. [10] and Nagano et al. [11,12] have published measurements 
of the fluorescence yield using electrons from the /3 decay of a ^^Sr radioactive 
source (average energy of about 0.85 MeV). Experiments at higher energy 
are performed using linear accelerators [10], [13,14,15]. Finally, measurements 
at low energy are being carried out with home made electron guns [16]. More 
information on the status of all these experiments can be found in [17]. On the 
other hand, the effect of uncertainties in the fluorescence yield on the shower 
reconstruction are being studied [18,19]. 



Electrons passing through the atmosphere deposit energy due to inelastic col- 
lisions with air molecules. Only a very small fraction of these processes (i.e. 
excitation of B^S+ state of N2' and C^n„ of N2) gives rise to the production of 
fluorescence light in our spectral range. However it is commonly assumed that 
the fluorescence intensity is proportional to deposited energy. This assumption 
has not yet been established neither theoretically nor experimentally, at least 
in a large energy range. 



In this paper, the air fluorescence yield is calculated using well known theo- 
retical properties and experimental data on the various cross sections involved 
in the fluorescence emission. Particular attention has been paid to the role of 
low energy secondary electrons ejected in molecular ionizations. As a result 
of our analysis recipes are given which allows the calculation of the absolute 
value of the total fluorescence yield (pressure and temperature dependent) as 
a function of electron energy in a wide interval ranging from the threshold to 
the Gev region. Comparison with available experimental data will be shown. 



2 Theoretical considerations 



In this section well-estabhshed properties on the molecular excitation (by elec- 
tron collision) and de-excitation are reminded, in particular, the relationship of 
the fluorescence yield with the excitation cross section and the transition prob- 
abilities (radiative and collisional). The contribution of secondary electrons to 
the fluorescence emission is evaluated. The relationship between deposited 
energy and fluorescence yield is discussed. All these properties will allow us 
to make theoretical predictions on the energy dependence of the fluorescence 
yield. 



2.1 Fluorescence yield and optical cross section 



Upon passage of an energetic electron through a molecular gas, the fluores- 
cence yield for a band e^y/ is deflned as the number of fluorescence photons 
emitted in the molecular transition v — v' per incident electron and unit path 
length, where v and v' represent the vibrational quantum numbers of the upper 
and lower electronic states of the transition respectively. At very low pressures 
and ignoring secondary processes, e„„i is given by 

£„„, = Na^B''''' = Na,,. , (1) 



where N is the number of molecules per unit volume, a„ is the excitation cross 
section for the v level, B'"'" is the branching ratio (i.e. the ratio between the 
partial A'^'" and total A^ radiative transition probabilities of the upper level) 

Avv' 4^""' 

Tjvv' r r /n\ 

EAvv' A'" '' 



and a^^r is the so-called optical cross section. As is well known [20], the ex- 
citation cross sections a^ are also approximately proportional to the Franck- 
Condon factors qx^v, deflned as the overlapping integrals between the vibra- 
tional wavefunctions of the lower and upper levels of the transition which in 
this case are the X ground state of N2 and the v excited state of N2 or N2' 
(depending on the band system) respectively. 



Therefore 



aoo qx-.o B^^ 



(3) 



The above relation allows the calculation of any optical cross section from data 
on a reference transition (e.g. (Too)- In particular, for the nitrogen molecule the 
parameters qx-*v and B^"" are available in the literature [21,22]. 



2.2 Quenching, temperature dependence and pollutants effect 



At high pressure, molecular de-excitation by collision with other molecules of 
the medium play an important role (collisional quenching). The corresponding 
transition probability AJ! is proportional to the collision frequency and thus to 
the gas pressure P (assuming fixed temperature). Defining P^ as the pressure 
for which the probability of collisional quenching equals that of radiative de- 
excitation Al{P'^) = Al 






K{P) = K^. (4) 



so the fiuorescence yield can be expressed in this case by 

^ ' A-, + Al{P) l + P/Pi 

The characteristic pressure is given by 



VnMkT 1 



where M is the molecular mass, k is the Boltzmann's constant, T is the gas 
temperature and a"„„ is the cross section for collisional de-excitation. 



Collisional quenching enlarges the total transition probability A" = A^ + A^ 
and therefore the lifetime of the population of excited molecules r^ = 1/A" is 
shortened as compared with the radiative one r^ = l/A"^ as 
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with r^ = 1/A^. Therefore lifetime depends on pressure as 
IIP 

-(1 + 4). 



T'^i^P) T^ PI 

Both r^ and P^ can be measured in a plot of reciprocal lifetime versus pressure 
(Stern- Volmer plot). This is a very well established technique in use since 



many years ago for the experimental determination of radiative lifetimes and 
quenching cross sections. 

For gas mixtures like air, the above formulation is valid taking into account 
that in equation (1) and (5) N is the density of fluorescence emitters (e.g. 
nitrogen amounts a 79% of the molecular density in air) and the quenching of 
excited molecules is due to collisions with the various components. Neglecting 
other components apart from nitrogen and oxygen the P^ parameter for air is 
written as 






where M„ and M^ are the molecular masses of nitrogen and oxygen respec- 
tively, o"„„ and ano are the cross sections for collisional quenching of the excited 
nitrogen with nitrogen and oxygen molecules respectively. The parameters /„ 
and fo are the molecular number fraction of nitrogen and oxygen in air. 

In some media (e.g. plastic scintillators), excited molecules may decay by other 
additional processes like intersystem crossing or internal conversion [23] and 
thus the lifetime is further shortened. However for nitrogen (air) gas such 
internal processes do not take place and therefore they will be neglected in 
this analysis. 

Equations (6) and (9) give the dependence of fluorescence yield with temper- 
ature. Firstly, P^ shows a vT dependence (i.e. collision frequency). On the 
other hand, the quenching cross sections (T„„ and ano are expected to depend 
on the kinetic energy of the colliders and therefore a further contribution to 
the temperature dependence could be expected. Unfortunately no data on this 
dependence for air in our temperature range is available in the literature and 
thus this effect, which is not expected to be relevant, will be neglected in our 
analysis. In fact measurements from the AIRFLY experiment are consistent 
with quenching cross sections independent on temperature [14]. 

As is well known oxygen is responsible for the large quenching of air fluores- 
cence as compared with the smaller auto-quenching of pure nitrogen (o"„o >> 
o"„„). Contribution of other air components (e.g. water vapor) or pollutants 
can be easily evaluated with equations (5) and the appropriate extension of 
(9) as far as the corresponding quenching cross sections are available. In this 
paper quenching by other air components will not be studied. 



2.3 Contribution of secondary excitations 



Apart from direct excitation by electron collision, the upper level can be popu- 
lated by de-excitation of higher energy electronic states. This cascading effect 
contributes to the fluorescence yield enlarging the apparent optical cross sec- 
tions, so experiments devoted to determine excitation cross sections need to 
make the appropriate corrections to the measured apparent values. While for 
our problem, this apparent cross section is the parameter of interest, no sig- 
nificant cascade effects have been reported and so this will be neglected in this 
analysis. 

Pure vibrational transitions (i.e. inside a given electronic state) induced by 
collisions with other molecules of the medium can redistribute somewhat the 
population of the upper levels. Very scarce data on this process are available 
in the literature. However the effect on the total fluorescence yield Y^w' ^w' 
(i.e. the total number of photons emitted per electron per meter) is expected 
to be negligible with a possible small effect on the spectrum of the fluorescence 
radiation. This effect will not be taken into account in our analysis. 

A very important secondary contribution will be that from low energy elec- 
trons ejected in ionizations. These excite molecular nitrogen increasing the 
observed fluorescence [24]. As aforementioned, for each primary incident elec- 
tron of kinetic energy i5, the number of direct excitation processes to the v 
level per unit length is N(7^{E). On the other hand, the number of generated 
secondary electrons per unit length is Naion{E) where aion{E) is the ioniza- 
tion cross section. For our calculations we will use the so-called gross ionization 
cross section which also includes processes in which more than one electron is 
generated in a primary interaction (even though contribution from those pro- 
cesses is small). If we name av{E, P) the fraction of these secondary electrons 
that excite the upper v level [ay^/ = a^B^^ the fraction that excites v — v' 
band emission), then an effective optical cross section al^l can be defined as 

<i' = (^vv'{E) + a,,iE, P)a,UE) , (10) 



and the total fluorescence yield given by equation (5) becomes 



The effect of these secondary electrons results the dominant mechanism for 
the production of 2P fluorescence at high pressure [24]. In addition, as will be 
shown in sub-section 3.1, a non- negligible contribution of the IN fluorescence 
is due to these secondary electrons. 



The a^' parameter can be calculated by means of a Monte Carlo simulation. 
For this work the simple algorithm described in [24] has been used with a 
somewhat improved treatment of the electrons escaping the interaction region. 
The resulting small dependence of a^^i on primary electron energy, pressure 
and characteristic size of interaction region i?, can be approximated (for IhPas 
and higher pressures) by 



avv'{E, P) = mmjsoln , Coln — | . 

si ei 



(12) 



This expression is purely empirical (a fit to the MC results). The value of the 
sq, Si, Co and ci parameters depend on the particular case. Results for 0-0 
bands of IN and 2P systems are shown in Table 1. 

The two terms inside the min function in equation (12) can be easily inter- 
preted as corresponding respectively to those conditions where electrons can 
escape from the interaction region (low P x R values) and those where elec- 
trons lose all their energy inside the interaction region (high P x R values). 
While in the first case fiuorescence by secondary electrons is limited by the 
traversed material, in the second one no (primary or secondary) electron es- 
capes from the interaction region and fiuorescence contribution is given by 
the primary incident energy. In both cases the (very slow) logarithmic de- 
pendence arises from the very small dependence of the secondary-electrons 
energy-distribution on the primary electron energy. For a given primary en- 
ergy, the saturation PR value (for which conditions change from one to the 
other behavior) can be interpreted as the range of secondary electrons in the 
gas and, as can be easily checked, this is almost independent of the studied 
process (IN or 2P fiuorescence). A paper describing more details of this model 
is in preparation [25]. 

Since a^^/ is proportional to the corresponding optical cross section, expression 
(3) leads to 



a,: 



aoo 






(13) 



which allows the calculation of the various a^y'^E, P) parameters from aoo- 





So 


si (hPaxcm) 


eo 


ei(eV) 


IN (0-0) 


8.67- 10-3 


1.92-10-2 


1.37-10-2 


73.7 


2P (0-0) 


2.00-10-3 


1.36-10-5 


3.22-10-3 


0.942 



Table 1.- Vaiues of the parameters in equation (12) for the most prominent bands 



of the IN and 2P molecular systems. 



Note that experiments determining quenching cross sections from the pressure 
dependence of fluorescence intensity may overestimate the P' values if they 
ignore the effect of secondary electrons by using equation (5) instead of the 
(11) right one. As mentioned above the P' parameter can be also measured 
in a plot of reciprocal lifetime versus pressure. This is a safe technique not 
affected by secondary electrons. 



2.4 Stopping cross section and deposited energy 



The stopping power for electrons S, i.e. the energy loss per unit length of 
traversed matter due to both excitation and ionization processes, can be ac- 
curately determined for energies above 1 keV by the well known Bethe-Bloch 
formula [26]. At lower energy, experimental results are available down to 25 
eV [27]. On the other hand, in general, S can be expressed [20] as 

S=^ = NY.En<Jn, (14) 



where the summation includes the energy En and cross section (T„ for all kine- 
matically accessible excited states n (ionization processes contributing there 
with a large family of continuum states). The cross section o"„ depends on the 
electron incident energy, and can be determined in electron scattering exper- 
iments [28]. A stopping cross section (XstiE), proportional to S, can be easily 
defined [20] by rewriting (14) as 

S = ^ = NRy Y.{EjRy)cXn = NRyast , (15) 



where Ry is the Rydberg constant (13.606 eV). This ast{E) cross section, 
proportional to the stopping power, will be useful for comparison purposes in 
the present analysis. 

The fluorescence efficiency ^^v' [11?12], [29] is defined as the fraction of de- 
posited energy which is emitted as photons of the v — v' band. Slow electrons 
deposit all the energy locally and thus the Bethe-Bloch formula gives the de- 
pendence of deposited energy on E. Therefore for low energy electrons 

™'^^ " ^^V^ " Rya^t l + P/Pl ■ ^ ' 
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Fig. 1. Cross sections (atomic units) involved in the emission of Euorescence as a 
function of energy for air at atmospheric pressure. The optical cross section ctqq^ for 
the IN system has been obtained from experimental data at low energy assuming a 
Born-Bethe behavior at high energy and including the density correction (see 3.1). 
The ionization cross section aion has been taken from [49]. Experimental results 
have been used for the 2P optical cross section a^Q (see 3.2). The agt cross section 
(proportional to the electron stopping power) has been obtained from [26] and [27] 
for energies above and below 1 keV respectively. Note that ast and aion follow a 
similar behavior although they are not exactly proportional. 

Note that the interpretation of the fluorescence efficiency given by some au- 
thors [1,11] is different from that of equation (16). These authors assume that 
all the electron energy loss is employed in the excitation of the molecules to 
one single upper level, a fraction is supposed to be emitted as fluorescence ra- 
diation and the remaining energy is lost by internal quenching. Such a model 
which is valid in some media (e.g. plastic scintillators) [23] is not applicable 
in our case. In fact, only a small fraction of the lost energy is used in the 
excitation of the level v of interest while the remaining energy gives rise to 
other direct excitations, not internal processes. The authors of ref. [11] argued 
recently along similar lines [29]. 

At high energy the stopping power also accounts for the production of high 
energy secondary electrons which do not deposit their energy locally (large 
range and bremsstrahlung emission). Unless the interaction region (e.g. gas 
cell in a lab experiment) is very large, a non-negligible fraction of the energy 
transferred to the medium by the primary electron is not deposited inside the 
volume and consequently a non-negligible fraction of the fluorescence radiation 
is not emitted in the interaction region. 



According to equation (11), several terms contribute to the total fluorescence 
yield. At low energy and low pressure the main contribution is direct excita- 
tion. In principle the energy dependence of agt and a^^i are not expected to be 
the same since o"„ excitation cross section only accounts for one of the many 
inelastic processes involved in the energy loss. On the other hand, even if o"„ 
followed the same behavior for all processes (including those leading to the 
generation of fluorescence light) e^' is not expected to be proportional to ast 
because of the different weighting values of En in equation (15) depending on 
each particular process. In Figure 1 the cross sections involved in fluorescence 
emission and energy deposition are shown as a function of energy for compar- 
ison. For instance, the different behavior of agt and 2P optical cross section is 
evident. 
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Fig. 2. Comparison on a Fano plot of the stopping cross section (Bethe-Bloch) and 
the ionization cross section (Born-Bethe) for air at IPa, 100 hPa and atmospheric 
pressure. Cross sections (in arbitrary units) have been normaUzed in order to get 
a common value at 1 MeV for all cases. Pure Born-Bethe function (e.g. ionization 
cross section at IPa ) is linear while the Bethe-Bloch stopping cross section clearly 
deviates from linearity. Density effect tends to saturate both functions at high energy 
although the effect is more pronounced for the ionization cross section. 

As shown in detail later, at high energy and high pressure the fluorescence yield 
is approximately proportional to the ionization cross section which follows the 
Born-Bethe relativistic law [20]. The Born-Bethe law for the electron-molecule 
collision can be expressed as 



a 



P'- 



■{\nC(3''-\n(l-/3^)-(3^} 



(17) 



10 



where (3 is the relativistic speed and, A and C are constants. At very high 
energies, a density correction has to be apphed to take into account the polar- 
ization of the molecules along the electron path [30,31,32]. For that purpose, 
we have included the density correction in the transversal contribution of the 
Born-Bethe expression [20], in the form proposed by [30]. The resulting ex- 
pression turns out to be 

a = ^{lnCp'-Hl-p')-P'-5F}, (18) 



where Sp is the usual density correction parameter included in the Bethe- 
Bloch formula for the stopping power [33] . Because of this density correction, 
the ionization cross section of the nitrogen molecule at high energy depends 
on pressure. 

For any pure Born-Bethe cross section (17), crf]"^ versus x = ln(/3^) — ln(l — 
jS"^) — /5^ is a straight line (the so-called Fano plot). In Figure 2 the ionization 
cross section has been represented on a Fano plot for several pressures together 
with the corresponding stopping cross section (Bethe-Bloch formula) [34]. As 
expected, at very low pressure, crjo„ follows a linear behavior while ast clearly 
deviates from linearity. At high pressure and very high energy both functions 
flattens as expected from the density effect but following a slightly different 
behavior. On the other hand Figure 2 shows that at low energy ast is not 
exactly proportional to cTion- 

In summary we can conclude that fluorescence yield is not expected to be 
exactly proportional to deposited energy since most of the fluorescence radi- 
ation is generated by low energy electrons (either primaries or secondaries) 
for which fluorescence emission (ctiat, o-2p) and stopping power {(Tion) follow 
different energy behavior. 



3 The fluorescence yield in pure nitrogen and dry air 



In this section experimental data on the excitation/de-excitation of the upper 
levels of the IN and 2P systems of nitrogen will be analyzed. This information 
will allow us the evaluation of the fluorescence yield. For all the calculations 
shown in this section a temperature of 300 K has been assumed. The density 
of fluorescence emitters (nitrogen molecules) A^ in dry air has been assumed 
to be 79% of the total molecular density. 
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3. 1 The First Negative system 

The IN system of Ng" has been studied by many authors. In particular many 
measurements of the optical cross section ctoo corresponding to the strongest 
band at 391.2 nm are available in the literature. At very low energies, absolute 
measurements have been reported by Borst and Zipf [35] (19 eV < E < 3.0 
keV), Stanton and St John[36] (19 eV < E < 400 eV), Aarts et al. [37] (97 eV 
< E <6 keV), Srivastava and Mirza [38] (70 eV < E < 4 keV), Holland [39] 
(90 ev < E < 2 keV), McConkey and Latimer [40] (26 eV < E < 300 eV) and 
McConkey et al. [41] (18 eV < E < 2 keV). In Figure 3 a significant set of 
the above measurements have been represented on a Fano plot. Also, relative 
measurements published in [42], normalized to an average value of the above 
authors in the overlapping region, have been included. 
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Fig. 3. Fano plot for the optical cross section croo (in atomic units) of the IN system. 
Experimental results (symbols) show a linear behavior at energies over about 150 
eV. Discrepancies in the slope are very likely due to calibration factors. An average 
linear fit (continuous line) has been used in this work for the calculation of the IN 
Euorescence yield. 

At higher energies some measurements are available. Data of O'Neil and 
Davidson in the range 3-60 keV [43] are shown in the figure. Cross sections 
inferred from previous measurements of these authors with 50 keV electrons 
at high pressure [44] are significantly smaller and they will not be used in 
the present analysis. As shown in the figure, the above measurements follow 
a linear behavior in a Fano plot. Note that the density correction is only rele- 
vant at very high energy. The slope of the various data sets, i.e. A parameter 
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in equation (17), have slightly different values while a reasonable agreement 
in the abscise in the origin is found, i.e. C parameter. This features can be 
interpreted as discrepancies in the absolute calibration factors of the various 
experiments, while their results agree in the shape of the o"oo(-E) function. 



An average linear fit of the available data to equation (17) is shown in Figure 
3. The values of the parameters were found to be A = 2.230 x 10~^ a^ and 
C = 9.826 X 10^. None of these experimental works reports any contribution of 
cascades to the measured optical cross section. As a further test of the absence 
of cascades, it has been checked [24] that the measured optical cross sections 
of this system fulfills equation (3). Therefore the above Born-Bethe function 
provides us with a reliable value of aoo in the whole energy range (see Figure 
1). Note that at very high energy the density correction has to be included. 
Optical cross section for other IN bands can be easily obtained using relation 
(3). 



While at very low pressures and low energies a^^/ values suffice for evaluat- 
ing fluorescence yields by means of equation (1), at those conditions we are 
interested here, quenching and secondary-electron contributions have to be 
considered by means of equation (11). Expression (12) and Table 1 show that 
the aoo parameter, accounting for the contribution of secondary electrons of 
the most intense band of the IN system (391.2 nm), ranges from (at low E 
and/or low PxR conditions) to 0.12 (at E = IMeV and PxR=5 atmxcm). In 
order to appraise the importance of these secondary contributions it must be 
noted that aoo — 0-075(rion{E) for IN system at large energies (see Figure 1). 
Introduction of that result in equation (11) reveals that secondary/direct con- 
tributions are in the q;oo/0.075 ratio and thus, contributions from secondary 
electrons can be very important at high pressures and high energies. 
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Fig. 4. Fluorescence yield of the most prominent band of the IN System (391.2 nm) 
for pure nitrogen at several pressures. Predictions of this work (continuous line) are 
compared with experimental data of Nagano et al. [12] (•) at the same pressures. 
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Fig. 5. Same as Figure 4 for dry air. In addition, experimental values of Hirsh et 
al. [45] at 20 hPa (*) and atmospheric pressure (o) are shown for comparison (see 
the inset for a more detailed comparison). 



The above results allow the calculation of eoo as a function of electron energy, 
pressure, temperature and characteristic size of the interaction region in a 
given experimental arrangement. As an example Figures 4 and 5 show eoo(E) 
for pure nitrogen and dry air respectively, in a wide energy range (eVs - GeVs) 
for several pressures (2, 20 and 1013 hPa) and a temperature of T = 300 K. 
Characteristics pressures of Pq = 1.7 and 1.3 hPa for pure nitrogen and dry 
air respectively have been assumed (see sub-section 3.3 for a discussion on 
available P' values). In these calculations a characteristic size oi R = 2.5 
cm has been assumed, which can be considered as representative for typical 
lab measurements. The fluorescence yield results of Hirsh et al.[45], for dry 
air, and Nagano et al. [12], for both nitrogen and dry air, are shown in the 
flgures for comparison. Direct measurements of [45] are in agreement with 
our predictions showing that our correction for secondary electrons accounts 
for these experimental results. On the contrary, our calculations are in clear 
disagreement with measurements of [12] for which a two line analysis had to 
be applied in order to subtract the contribution of other band not separated 
by the spectroscopy fllter. As shown below, the contribution of the IN system 
to the total air fluorescence yield (at high pressure) is in the range of 15% and 
therefore the above discrepancy is not so relevant for the total fluorescence 
yield. 
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3.2 The Second Positive system 



Fluorescence of the 2P system in both pure nitrogen and air has been exten- 
sively studied by many authors. Experiments aiming at measuring the electron 
excitation cross section of the upper level C^IIu carried out on pure nitrogen 
at low electron energy {E < 1 keV) and low pressure (a few mTorr or lower) 
show a fast E'"^ dependence of the fluorescence light [46,47,48] much faster 
than that of the IN system. This result is expected, taking into account the 
optically-forbidden nature of the 2P transition. In Figure 1 the optical cross 
section for the — band of the 2P system from the abovementioned mea- 
surements is shown. 
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Fig. 6. Fluorescence yield of the most prominent band of the 2P System (337.0 nm) 
for pure nitrogen at several pressures. Predictions of this work (continuous line) are 
compared with experimental data of Nagano et al. [12]. 

Since at 1 keV energies, 2P fluorescence observed at low pressure is much 
weaker than that of the IN system, a naive extrapolation would predict that 
the 2P fluorescence should be completely negligible compared to IN one at 30 
keV energies and above. On the contrary experiments [10,11,12], [44] carried 
out at high energies (up to 1 GeV) on air at high pressure (ranging from a 
few hPa up to 1 atm) shows that the 2P fluorescence even dominates over 
the IN system. This feature can be explained taking into account the effect of 
the secondary electrons ejected in ionization processes [24]. The contribution 
of these secondary electrons have been evaluated from a MC calculation as 
explained above and the values of the corresponding parameters of expression 
(12) for the most intense band of this system (337.0 nm) are shown in Table 
1. 

Characteristic pressures of Pq = 77 and 13.1 hPa for pure nitrogen and dry 
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Fig. 7. Same as Figure 6 for dry air. In addition experimental results of Ueno [13] 
at atmospheric pressure (o) are shown. 

air respectively have been assumed (see below for a discussion on available P^ 
values) . Figures 6 and 7 show the predicted fluorescence yield for the 0-0 band 
of the 2P system as a function of electron energy from threshold (eVs) up to the 
GeV region at several pressures in pure nitrogen and dry air respectively. An 
interaction region of 2.5 cm radius has also been assumed for this calculation. 
Experimental data for the 337.0 nm band at 0.85 MeV [12], 1 MeV and near 
1 GeV [13] are shown for comparison. 



3.3 Total fluorescence yield 



The total fluorescence yield, i.e. the sum of the contributions from all the 
molecular bands of both systems IN and 2P, 









(19) 



can be easily calculated from (11) provided that the involved molecular pa- 
rameters are available for all transitions. Taking into account (3) and (13) 
the contribution of each band system to the total fluorescence yield can be 
expressed as a function of a reference transition (e.g. 0-0) as 
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where (2) has been used. All the parameters depend on the system, IN or 2P. 
The terms inside the v-v' sum represent the relative intensities of the bands 
for the corresponding system and at high pressures (P ^ P^) they can be 
approximated by 



e..' ~ gx^.5™'P: . 



(21) 



In general, equation (21) is able to predict experimental results on relative 
intensities although a detailed comparison is beyond the scope of this paper. 



v/ v' 





1 


2 


3 


5 


5 


qx^v 


PU (hPa) 


P„t (hPa) 




391.2 


427.5 


470.6 


522.5 


586.1 


665.9 











0.627 


0.204 


0.043 


0.007 


0.001 


0.000 


0.883 


1.7 


1.3 




358.0 


388.2 


423.4 


464.9 


514.6 


575.0 








1 


0.041 


0.029 


0.030 


0.011 


0.003 


0.001 


0.114 


- 


- 




330.5 


356.1 


385.5 


419.7 


459.7 


570.4 








2 


0.000 


0.001 


0.000 


0.000 


0.000 


0.000 


0.002 


- 


- 



Table 2.- Molecular parameters for the 1 N system. Each box of the vv' table 
shows the wavelength (nm) of the v-v' transition (upper number) and the product 
B'"'" ■ qx~iv (lower number). The horizontal sum of these products is equal to the 
qx^v Franck-Condon factor for the molecular excitation. Last two rows show the 
values of the characteristic pressures for the quenching of the upper v level in pure 
nitrogen and air (see text for more details). 



v/ v- 





1 


2 


3 


4 


5 


6 


7 


8 


qx^v 


Kitr (hPa) 


P^.^ (hPa) 




337.0 


357.6 


380.4 


405.8 


434.3 


466.5 


503.2 


545.2 


593.8 











0.265 


0.175 


0.072 


0.022 


0.006 


0.001 


().()()() 


0.000 


0.000 


0.545 


77 


13.1 




315.8 


333.8 


353.6 


375.4 


399.7 


426.8 


457.3 


491.7 


530.9 








1 


0.138 


0.007 


0.064 


0.057 


0.028 


0.010 


0.003 


0.001 


0.000 


0.308 


36 


11.2 




297.6 


313.5 


330.9 


349.9 


370.9 


394.2 


420.0 


448.9 


481.3 








2 


0.016 


0.041 


0.003 


0.007 


0.016 


0.013 


0.006 


0.002 


0.001 


0.106 


23 


9.1 




281.8 


296.2 


311.5 


328.4 


346.8 


367.1 


389.4 


414.0 


441.5 








3 


0.001 


0.009 


0.007 


0.003 


0.000 


0.003 


0.004 


0.002 


0.001 


0.030 


22 


7.9 




268.4 


281.2 


295.2 


310.2 


326.6 


344.5 


364.1 


385.6 


409.3 








4 


0.000 


0.000 


0.003 


0.001 


0.001 


().()()() 


().()()() 


0.001 


0.001 


0.008 


21 


6.7 



Table 3.- Same as Table 2 for the 2 P system. 

Values for the involved parameters are shown in Tables 2 and 3 for IN and 
2P systems respectively. Accurate Franck-Condon factors qx-^v and branch- 
ing ratios B™ have been taken from [21]. In regard with the characteristic 
pressures P^, many measurements are available in the literature. 
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For the IN system [1], [12] and [50] provide experimental results for the v=0 
level for both pure nitrogen and dry air. A reasonable agreement between [1] 
and [50] is found and an average value of both results, shown in Table 2, has 
been used in our calculations. Since no results for f > seems to be available, 
the same Pq value has been used for all v levels. The bands of this system 
are strongly attenuated by quenching and therefore their contribution at high 
pressure is very small, consequently, even a severe error in P^ {v > 0) would 
not have a significant effect on the calculated total fluorescence yield. 

On the other hand, quenching of the 2P system by nitrogen (auto-quenching) 

has been extensively studied by many authors [1], [12], [50], [53,54,55,56,57,58,59,60,61,62,63,64]. 

Unfortunately large discrepancies between the various measurements are found. 

Averaged values from [59], [60] and [61], in good agreement with a significant 

fraction of previous measurements, have been used in this work (see Table 3). 

The total fluorescence yield for the bands contributing to the usual experimen- 
tal spectroscopic interval (300 - 406 nm) has been evaluated as a function of 
electron energy and pressure. Note that expression (20) had to be corrected to 
account for the missed transitions (those outside the above wavelength limits). 
A plot of e{E) for pure nitrogen in the interval ranging from threshold to the 
GeV region at several pressures ranging from 1 hPa to atmospheric pressure 
is shown in Figure 8. Several interesting features of e{E, P) can be observed 
in this figure. At high pressure, expression (20) becomes 

system _ ^ ^00 (-^; "j y^ p/ fcycy\ 

showing that for P ^ P^ the fluorescence yield is basically independent on 
pressure. The only remaining dependence is the small one from the a parame- 
ter and (even smaller) that from the density correction of (Tion at high energy. 
That is, the increase of fluorescence with pressure is nearly canceled out by 
the corresponding increased quenching. Note that the higher is P^, the weaker 
is the pressure dependence of fluorescence yield. 

The o"oo parameter can be easily evaluated. For energies about 100 eV and 
over, the optical cross sections of the IN system are proportional to the total 
ionization cross section, ctIq^E) = Xoo'^ion{E) and therefore 

al^'^ff{E, P) = (xoo + al^)cr^on • (23) 

In regard with the 2P system the contribution of direct excitation at high 
energy and high pressure {E > lO^eV and P x R> 3hPa x cm) is negligible 
because under these conditions a^^i <^ Qw'O'ion and thus 

ali;'^ff{E,P)^al^a,on. (24) 
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Fig. 8. Total Euorescence yield (300 - 406 nni) versus energy for pure nitrogen at 
pressures (bottom - up) 1, 2, 10, 20, 100 and 1013 hPa for a characteristic size of 
the interaction region of R = 2.5 cm. 

The a(E,P) parameter reaches saturation (i.e. it becomes independent of E 
for increasing energies) at an energy which grows with pressure (12). The data 
in Table 1 shows that, at atmospheric pressure, this saturation energy results 
of about 1.3x10^ eV for both systems. As a consequence for higher energies 
e(E) is proportional to aion{E). 

The results of our calculations have been compared with available experimen- 
tal data. Figure 9 shows the total fluorescence yield (300 - 406 nm) at high 
energy and atmospheric pressure for pure nitrogen together with the absolute 
measurements of Nagano et al. [12] at 0.850 MeV and those of the FLASH 
collaboration [15] at 28.5 GeV. In addition, relative measurements of AIRFLY 
[14] in the interval 50 - 500 MeV have been represented. As shown in the fig- 
ure, experimental data for pure nitrogen are in good agreement (within the 
experimental errors) with our predictions. 

Unfortunately, quenching measurements in air for the 2P system are scarce, 
showing also large discrepancies. Bunner [1] gives P^ values of 20 and 8.7 hPa 
for f = and 1 respectively. Nagano et al [12] measure P' values from the 
fitting of experimental £{P) results. Average values of about 18, 25 and 8 
hPa are found for f = 0, 1 and 2 respectively. As already mentioned these 
P^ values could be overestimated by the contribution of secondary electrons. 
Finally, Pancheshnyi et al. [60] have published measurements of the Nitrogen 
2P quenching by O2 molecules in good agreement with the results of [64] and 
[55]. These measurements give the air P^ values shown in Table 3 which have 
been used in our calculations. 

Figure 10 shows plots of £{E) for dry air in a wide energy interval for sev- 
eral pressures ranging from 1 hPa to atmospheric pressure. In Figure 11 these 
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Fig. 9. Total fluorescence yield (300 - 406 nni) at high energy for pure nitrogen. Pre- 
dictions of this work (continuous line) are compared with absolute measurements of 
Nagano et al. [12] (•) and the FLASH collaboration [15] (k). Relative measurements 
of AIRFLY (o) normalized to our predictions at Iff eV are also shown. Dashed line 
represents the Bethe-Bloch function normalized at 1 MeV. 

results are compared with available measurements at high energy. Our calcu- 
lations are in good agreement with that of Nagano et al. [12] at 0.850 MeV 
but they predict a fluorescence yield about 20% larger than the measurements 
of Kakimoto et al. [10] (at several energies in the range 1 MeV - IGeV) and 
those of the FLASH collaboration [15] at 28.5 GeV. 

The calculations presented in this work relies in several molecular parameters 
and in a model for evaluating the contribution of secondary electrons. Branch- 
ing ratios, Franck-Condon factors and the validity of the Franck-Condon prin- 
ciple are reliable and thus their contribution to the uncertainty in e is ex- 
pected to be very small. On the other hand, at high pressure, fluorescence 
yields are proportional to P^ values for which large discrepancies are found 
in the available measurements. The uncertainty in these parameters are very 
likely one of the main error sources in our calculations. On the other hand, our 
model [24,25] uses several molecular parameters which can also contribute to 
a systematic error in the evaluation of the fluorescence light from secondary 
electrons, in particular, the absolute value of the total excitation cross section 
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Fig. 10. Total fluorescence yield (300 - 406 nm) for dry air against electron energy 
at pressures (bottom-up) 1, 2, 10, 20, 100, 1013 fiPa for a cliaracteristic size of the 
interaction region of R = 2.5 cm. 

of N2. This uncertainty leads to a possible (small) constant factor error. In 
addition, the energy spectrum of secondary electrons may affect the energy 
dependence of the fluorescence yield. In this work the spectrum used in [24] 
has been assumed. On the contrary, the geometry of the interaction volume 
(different for the various experiments) should not contribute significantly to 
our errors since «„„' grows very slowly with R (^ InR). In regard with the 
energy dependence, our calculations rely in the validity of the well established 
Born-Bethe approximation at high energy. 

In summary, at high pressure the predictions of this work could be affected by 
a small error (constant factor) mainly due to the uncertainty in the quench- 
ing efficiency together with a possible error in the energy dependence due to 
inaccuracies in the spectrum of secondary electrons. In Figure 11 our fluo- 
rescence yield reduced by a constant factor of 18% has been plotted showing 
a reasonable agreement with all experimental results. Note that the ratio of 
fluorescence yields 5««»'/5™*''05en ^^ experiments measuring both gases (air and 
nitrogen) at high pressure are expected to be equal to the ratio P'airl PLitrogen 
since experimental errors (calibrations, geometrical factors, etc.) are the same 
for both gases. Our predictions for nitrogen agree with the absolute measure- 
ments of [12] and [15] while a fluorescence efficiency about 20% larger than [10] 
and [15] is predicted for air. This lead us to suspect that the ratio PUrl ^nitrogen 
assumed in our calculations has been overestimated by about a 20%. 

In Figures 9 and 11 a plot of the stopping power law for electrons (Bethe- 
Bloch) has been represented showing a noticeable deviation with respect to 
s{E) as predicted, since at large E values a significant fraction of the electron 
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Fig. 11. Total Euorescence yield (300 - 406 nni) at high energy for dry air. Contin- 
uous line represents the prediction of this work using the parameters of Tables 2 
and 3. The comparison with experimental results of Nagano et al. [12] (•), Kaki- 
moto et al. [10] (o) and the FLASH collaboration [15] (k) suggests that P^ values 
could be overestimated (see text for more details). Results of this work reduced by 
18% (dash-dotted line) are in good agreement with all measurements. Dotted line 
represents the stopping power and the dashed line is the deposited energy from [15] 
both normalized to e at 1 MeV. 



energy loss is not deposited inside the fluorescence cell. At low energy the 
fluorescence yield also deviates from the stopping cross section as shown in 
Figure 12. As already mentioned most of the fluorescence light is generated 
by low energy secondary electrons for which deposited energy is proportional 
to the stopping cross section. Since at low E, e/agt is energy dependent, the 
e to deposited energy ratio is expected to depend on E at high energy also. 
In Figure 11 the deposited energy versus E reported by [15] for the FLASH 
experiment has been also represented (normalized to our fluorescence yield 
at 1 MeV). Our calculations predict that at large energies a growing fraction 
of deposited energy is converted in fluorescence radiation. According to this 
flgure, the ratio e to deposited energy grows about 19% in the interval 1 MeV 
- 1 GeV. 
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Fig. 12. Ratio of Euorescence yield over stopping cross section (Bethe-Bloch) for 
air (arbitrary units) at 1 hPa (dotted), 20 hPa (dashed), 200 hPa (dot-dashed) and 
1013 hPa (continuous). The ratio has been normahzed to 1 at IMeV. 

4 Conclusions 



Theoretical predictions on the fluorescence generation of nitrogen molecules 
excited by electrons have been presented. A general procedure for the cal- 
culation of fluorescence yields is shown which can be applied in a very wide 
energy interval ranging from threshold up to the GeV region and for any envi- 
ronmental condition (pressure, temperature, contaminants, etc.) as far as the 
involved molecular parameters are available. As an example, plots of fluores- 
cence yields (as expected from typical laboratory experiments) versus energy 
in a wide pressure range are shown. 

The contribution of direct excitation of IN system is calculated using low 
energy measurements of the optical cross section with extrapolation to high 
energies as given by the Born-Bethe law. Direct excitation of 2P system is 
negligible. Particular attention has been paid to the the contribution of sec- 
ondary electrons ejected in ionization processes. A simple model is used to 
evaluate the fluorescence from secondary electrons which turns out to be non- 
negligible for the excitation of the IN system and dominant for the generation 
of the 2P fluorescence. At high energy and high pressure fluorescence yield 
results proportional to the total ionization cross section which also follows a 
Born-Bethe energy law. Quenching of fluorescence by N2 and O2 molecules 
is taken into account using available data in the literature. At high pressure, 
fluorescence yield is proportional to the P^ parameters. Unfortunately values 
of these parameters reported in the literature show large disagreements. 

Comparison of our predictions with available measurements at high energy 
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{E > IMeV) and high pressure shows very good agreement for pure nitrogen 
while some discrepancies of about 20% are found for dry air. These discrep- 
ancies are very hkely due to uncertainties in the quenching cross sections. 

The relationship between fluorescence yield, stopping power and energy de- 
posited by the electrons in the medium has been discussed. In a typical lab- 
oratory experiment, at very high energy, many secondary electrons reach the 
wall of the cell before emitting all the fluorescence light and therefore stopping 
power grows with energy faster than fluorescence yield. On the other hand, 
most of the fluorescence light is generated by low energy secondary electrons 
for which the optical cross section is not proportional to the stopping power 
and therefore, in principle, proportionality between fluorescence yield and de- 
posited energy is not assured. A comparison of our fluorescence yield results 
with deposited energy in the FLASH experiment shows a deviation from pro- 
portionality of about 20% in the interval 1 MeV - 10 GeV. Note that the 
uncertainty of this result depends on possible inaccuracies in the energy spec- 
trum of secondary electrons assumed in our model [24]. 

More details on the model used in this work for the evaluation of the fluo- 
rescence contribution of secondary electrons will be published soon. On the 
other hand a detailed study on the relationship between deposited energy and 
fluorescence yield following our model is underway. 
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